STM imaging of a bound state along a step on the surface of 
the topological insulator Bi2Te3 
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Detailed study of the LDOS associated with the surface-state-band near a step-edge of the strong 
topological- insulator Bi2Te3 , reveal a one-dimensional bound state that runs parallel to the step- 
edge and is bound to it at some characteristic distance. This bound state is clearly observed in 
the bulk gap region, while it becomes entangled with the oscillations of the warped surface band at 
high energy [l[, and with the valence band states near the Dirac point. We obtain excellent fits to 
theoretical predictions 01 that properly incorporate the three-dimensional nature of the problem to 
the surface state [3]. Fitting the data at different energies, we can recalculate the LDOS originating 
from the Dirac band without the contribution of the bulk bands or incoherent tunneling effects. 
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612X63 and 612863 hav6 b66n argu6d r6C6ntly to b6 
thr66-dim6nsional (3D) topological insulators (TI) 0-01' 
exhibiting a bulk gap and a single, non-degenerate Dirac 
fermion surface band topologically protected by time- 
reversal symmetry Subsequent angle resolved pho- 
toemission spectroscopy (ARPES) and by scanning 

tunneling microscopy (STM) and spectroscopy (STS) 
[H, [13, El experiments confirmed that prediction. Fo- 
cusing on 6i2Te3, ARPES revealed that with appropri- 
ate hole-doping, the Fermi level could be tuned to inter- 
sect only the surface states, indicating fully gapped bulk 
states as is expected from a three-dimensional TI [1, Q . 
Complementing ARPES, STM and STS studies empha- 
sizing quasi particle scattering interference from impu- 
rities and macroscopic defects, demonstrated the "pro- 
tected" nature of the surface band [H, [13, HI- 

While protection of the surface state is guaranteed 
within a simple Dirac band, warping effects can pro- 
duce special nesting wavevectors that allow for electronic- 
waves within the surface state band (SS6) These 
effects have been observed experimentally in the upper 
part of the SS6 by Alpichshev et al. In that pa- 

per we reported STS studies on high-quality 6i2Te3 crys- 
tals which exhibit oscillations of LDOS near a step-edge. 
While within the main part of the SS6 oscillations were 
shown to be strongly damped, supporting the hypothe- 
sis of topological protection, at higher energies, as the 
SSB becomes concave (as is observed by complementary 
ARPES data), oscillations appear which disperse with a 
particular wave- vector corresponding to the allowed spin 
states of the warped constant-energy contour [5| . 

In this paper we present a more detailed study of the 
LDOS associated with the SS6 near a step-edge, and 
report on the observation of what is most probably a 
one-dimcnsional (ID) state bound onto it. This bound 



state is clearly observed in the bulk gap region, while 
it becomes entangled with the oscillations of the warped 
SSB at high energy and with the valence band states 
near near the Dirac point. We successfully fit the data 
to theoretical predictions that properly incorporate 
the three-dimensional nature of the problem (outlined in 
[3) to the surface state . Fitting the data at different en- 
ergies, we can recalculate the LDOS originating from the 
Dirac band without the contribution of the bulk bands 
or incoherent tunneling effects. The excellent agreement 
between the experimental data and simple theory is an- 
other testimony to the simple, yet robust nature of the 
surface state in this model system of TI. 

For the present study we used Cd doped single crys- 
tals of Bi2Te3 identical to those used in our earlier study 
of the warped surface band [1]. Nominal doping levels 
of up to 1% for Cd were incorporated to compensate n- 
type doping from vacancy and anti-site defects that are 
common in the 6i2Te3 system. Actual doping was deter- 
mined separately using chemical and Hall-effect methods 
and were shown by ARPES Q to be in excellent agree- 
ment with the relative position of the Dirac point with 
respect to the Fermi energy. For the data described in 
this paper the Dirac point was located near ^300 meV. 

Samples were cleaved in vacuum of better than 5 x 
j^Q-io q^Qj-p^ a,nd quickly lowered to the ~9 K section 
of the microscope, where cryo-pumping ensures that the 
surface remains free from adsorbates for weeks or longer. 
Topography scans were taken at several bias voltages and 
setpoint currents (usually 200mV and lOOpA). A thor- 
ough discussion of the important features of the data that 
was collected from many topography and spectroscopy 
scans can be found in [1], here we concentrate on scans 
in the vicinity of a step-edge defect that was obtained in 
the process of cleaving the crystal. The measured thick- 
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ness of the step is ~30.5 A (in excellent agreement with 
the thickness of one unit cell [l^ ) , and we concentrate on 
the region away from the step as is shown in Fig. [TJ We 
mark the beginning of the step which runs along the [100] 
direction with x — (this choice is arbitrary to within a 
few A because of the roughness of the step) . 








FIG. 1: (color online) bottom: Topography of the step. The 
total height of the step ~30 A corresponding to one unit cell; 
top: DOS at E — —160 meV in the plateau region away 
from the step showing a bound state about 15 A away from 
the step. We also show in a cartoon the tip-step interaction, 
showing that pure surface LDOS can be taken on the plateau 
away from the step (A), but due to the strong decrease in 
height to the left of the step, rounding of the topography 
must occur (B). 

Above the topography, we also show in Fig. [T] the 
LDOS of that region for a bias voltage -160 meV (in 
the middle of the SSB) which lies in the gap between the 
bulk conduction band (BCB) and the bulk valence band 
(BVB). While the LDOS is rather flat and featureless in 
most of the scanned pronounced LDOS-peak is 

observed next to the step with constant strength along 
the step. We identify this peak as a bound state and 
study it further as a function of energy. 

As described in ref. [l|, at high enough energies, away 
from the Dirac point and in the region of warped surface 
band contour we observe oscillations of the LDOS that 
can be fit with the expression: p{x) oc A{E)sm{2kx + 
(j))/x (amplitude A, and a scattering potential phase-shift 
(/)). In this region LDOS demonstrates only simple and 
naturally expected deviation from the simple asymptotic 
expression around the origin. 

As the energy is lowered, and we enter the region 
in which oscillations are strongly damped (below ~ 
— lOOmeV), a new peak seems to emerge at a distance 
^ I5A from the step. For example, we show in Fig. [2] 
raw LDOS data as a function of distance from the step 
(zero marks the position of the step) in both direction, 
for bias energy of —190 meV, inside the bulk gap. For ref- 
erence we also show in the same figure a cross section of 
the topography of the step. Note that due to the height 



of the step ( 30 A) , when the tip drops down to keep the 
distance to the surface constant, it "probes" the edge of 
the step which results in a rounded region that extends 
to — 50A. We estimate that the actual step drops down 
to the lower region rather steeply around zero. To the 
right of the step the peak at ^15 A is the only feature 
that exists. The solid line shows a fit to our model [3], 
as will be discussed below. We note that we can impose 
the mirror-image of the fit to the "behind the step" re- 
gion. The excellent agreement that is observed beyond 
the region in which the tip did not clear the step is a 
direct proof of the two-dimensional nature of the sur- 
face band that beyond the exponentially decaying effect 
of the bound state, it "does not know" about the defect 
as if there is one, uninterrupted surface that wraps the 
crystal. 
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FIG. 2: (color online) LDOS as a function of position from 
the step at _E = —190 meV (in the region of the SSB that is in 
the bulk gap) shown together with a fit to Eqn.[4](solid line). 
Also shown is a cross-sectional profile of the step topography. 
In the shaded region the tip interacts with the side of the step 
and its lower part (see Fig. [1]), thus cannot be used for a fit. 

Fig. |3] shows the raw data of LDOS near the step with 
similar fits to Fig. [S] The result at —115 meV seems 
to be a smooth progression from higher energies, if we 
subtract the oscillatory part, where by —150 meV the 
peak stops changing shape and it can be fitted with our 
model with no residual. 

Despite the robustness of the experimentally observed 
peak, neither a clear prediction, nor possible explanation 
exists to date to account for this phenomenon. We at- 
tribute this deficiency to the pure-2D theoretical frame- 
work which all of the previous studies used. The results 
of these calculations (e.g. [11]) show that any features 
obtained are necessarily dispersive with energy while the 
peak in LDOS we report on is clearly not, which unam- 
biguously implies some details are missing. 

While the initial Hamiltonian derived from the band 
structure of Bi2Te3 did include all bands pi], all the calcu- 
lations of impurity scattering or scattering from macro- 
scopic defects used the truncated two-dimensional k x (t 
effective surface Hamiltonian first introduced by Zhang 
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FIG. 3: (color online) LDOS oscillations as a function of the 
distance from the step for several energies in the range of the 
SSB that is inside the bulk gap 1]. solid lines are fits to 
Eqn. 3] (see text). 



et al. 0. At most, higher order surface terms were intro- 
duced to account for the hexagonal warping effects [12 1. 
Here we claim that by using a pure 2D Hamiltonian, one 
inherently assumes that and ky are "good" quantum 
numbers, which may be a wrong starting point in the 
case of a strong perturbation such as a step-edge. When 
the local curvature of the surface of a topological insu- 
lator is much less than the characteristic length-scale of 
the problem ^ vp/^buik {^buik is the bulk gap), one 
cannot assume that the surface state will be just a so- 
lution of the conventional "k x cr" surface Hamiltonian, 
wrapping the new curved surface. Instead, effects due 
to the gradients near the corrugated surface induce bulk 
interference which become relevant. Therefore, a coher- 
ent treatment of such surface defect should start with the 
appropriate 3D Hamiltonian such as the one derived by 
Zhang et al. [3] (we choose the units such that fi — 1): 



nijir) = {M(r)l2x2®r^ 

+ v(k.j,crx ®Tx + ky(jy (g) Tx) + Vz{k^az ® t^) 

+ 0{k^)]ij{v) ^ E^j{v) (1) 

where Ui and ti are Pauli matrices in spin and orbital 
space respectively with i — x]y; z, and I2x2 is the 2x2 
identity matrix, v, Vz, and M(r) are material-specific 
parameters. In particular, the "mass term" M(r) rep- 
resents the asymmetry between the inside and outside 
the bulk of the TI material and has units of energy in 
our notations. Asymptotically M(r) = — Mi„ < inside, 
and M(r) = Mout > outside, while near the step M(r) 



varies rapidly to match the asymptotic behavior. It is 
easy to see that from the topology point of view, it is 
enough to set M(r) = +AI in vacuum and M{r) = —Al 
in the TI bulk in order to create a boundary that sep- 
arates the two distinct topological states . Further, to 
localize the effect of the mass twist to the TI surface the 
Eq. [1] can be squired. The result is a new Schrodinger- 
like equation (for convenience, the z-axis was re-scaled 
such that z = z{v/vz)) for which it can be demonstrated 
that the LDOS obtained from it is directly related to the 
LDOS of the original unsquared Hamiltonian. 



- v{VM{r) ■a)(g,Ty }i/)(r) = E^i;{r). (2) 

However, despite the fact that the effect of the step is now 
separated into the gradient term, Eq.[2]is still difficult to 
solve due to the unique spin structure associated with 
the eigenstates. However, one may argue that since at 
least one of the edges of the step (the one perpendicular 
to the cleavage surface) is very rough, it may serve as a 
termination line for the spin state. 

Close examination of the geometry of the step as shown 
in Fig. [1] indicates a bound surface state that is ^ 15 A 
away from the x = axis, and spreads along the y-axis. 
The step, which is a line defect in the x-y plane, is created 
by an intersection of a surface perpendicular to the z-axis 
with a surface perpendicular to the x-axis. In such a 
formulation ky continues to be a good quantum number 
in the plane. Therefore, setting ky = 0, we make an 
approximation in which we assume that we can ignore 
the spin structure on the rough surface of the step, hence 
obtaining the equation: 



[IVr -i;^V^)|V'| - w|VM| 







(3) 



Solving to the first approximation this equation without 
the spin variables as discussed above yields: 



MM 

1^0 I 



-x/A 



Koix'/X) 



dx' (4) 



where It/jqI is the asymptotic amplitude of the wave- 
function far away from the step, Kq{x) is the zero*''- 
order modified Bessel function of the second kind, and 
A = v/M ~ fiVF / ^buik — 30A is the characteristic 
length-scale for variations in the x-direction. 

An examination of the above solution shows that 
■0(a;)|/|'!/'o| has a single maximum {ipmax) near the step 
and then decays exponentially far away from the step. 
The corresponding LDOS ratio, a{E) = Pmax/po ~ 



/iV'ol , is found to be approximately 1.2, and 



independent of energy for the wavefunction solution in 
Eqn. m Obviously, the total density of states will be built 
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of a combination of wavefunctions that take into account 
the actual roughness of the surface and the correct spin 
configuration. While such a calculation is beyond the 
scope of the present paper, we expect that the shape of 
the correct solution solution, i.e. a density of states that 
peaks close to the step and decay asymptotically far from 
the step, will be of the same nature [ij. Apart from par- 
ticular details, within the presented approach the struc- 
ture of the peak in LDOS is determined by the bulk prop- 
erties of the samples, i.e. on the energy scales of the 
gap size implying that in the low energy regime (close to 
the Dirac point) the LDOS profile should be independent 
of energy. This is certainly the case for the position of 
LDOS maximum and the shape of the profile up to a ver- 
tical shift. However it also means that the relative height 
Pmax/ Pf) of the LDOS must be independent of energy too. 
Fig. |4] shows a typical LDOS spectrum away from the 
step, together with the ratio of peak height to asymptotic 
level of the LDOS as determined from the experiment as 
a function of bias voltage. What it shows though is that 
in the exposed region of the Dirac band, the ratio a[E) 
is clearly energy dependent. We assume that this depen- 
dence is due to an incoherent density of states Pinc{E), 
such that ptotai{E,x) = pcoh{E,x) + p.ndE). Assuming 
that the coherent part of LDOS pcoh{E, oo) behaves as 
explained above we can obtain from here an expression 
for it in terms of measurable quantities only: 

p,ok{E, oo) = pUE, oo) • f^JEA^)_l.\ (5) 

Here amax is the relative hight of the peak in the case 
the was no incoherent contibution to the LDOS. Since 
the latter can only decrease the height of peak, for speci- 
ficity we take amax = 1-35 here, which is the maximum 
value of aexp{E)- The actual value for amax should be 
given by the complete theory. The resulting coherent 
LDOS obtained this way is plotted in FigH] as a func- 
tion of energy. It is seen immediately that the resulting 
line is strikingly straight and points exactly to the po- 
sition where the Dirac point is expected to be located 
based on electron wave analysis presented in as well 
as ARPES data on these samples To justify that 
we extrapolated the validity region of EqlS] quite far way 
from the Dirac point although it is formally applicable 
only in the vicinity of DP, we point out that first of all, 
the actual criterion was (E/M)^ ~ (E/A)"^ « 1 which 
is still roughly met when {E — Ep) w —lAOmeV; and 
second it is seen from Fig. [3] that for these energies the 
Eqn. m is still providing a good fit, suggesting that the 
physics involved in obtaining Eqn. |4]is still valid. 

In conclusion we report on the observation of an ac- 
cumulation of DOS near an atomic step on the surface 
of Bi2Te3 identified as a bound state, whose structure is 
found to be in nice quantitative agreement with the the- 
ory developed in ref |2| . We also introduce a technique to 
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FIG. 4: (color online) LDOS spectra as a function of energy 
(left axis) away from the step (a), at the position of the max- 
imum of the peak (b), the ratio between the peak height and 
the asymptotic LDOS (c) , and the extracted coherent part of 
LDOS (d) which is proportional to the difference between the 
blue and green curves (see text). Also shown is the position 
of the Dirac point for this crystal inferred from the analysis 
of LDOS oscillations and ARPES. The region between the 
vertical dotted lines is within the bulk gap of the crystal. 



extract the coherent part of the DOS, and demonstrated 
that unlike the total measured DOS, its shape is in ac- 
cordance with the expected profile of DOS of a 2D Dirac 
band. 
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